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We examine the role of the condensing agent in the formation of polyelectrolyte bundles, via grand- 
canonical Monte Carlo simulations. Following recent experiments we use linear, rigid divalent cations 
of various lengths to induce condensation. Our results clarify and explain the experimental results 
for short cations. For longer cations we observe novel condensation behavior owing to alignment 
of the cations. We also study the role of the polyelectrolyte surface charge density, and find a 
nonmonotonic variation in bundle stability. This nonmonotonicity captures two trends that have 
been observed in separate experiments. 
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Under conditions commonly found in biological sys- 
tems, like-charged polyelectrolytes in solution can form 
compact aggregates. This phenomenon, referred to 
as polyelectrolyte condensation, occurs for a variety of 
biopolymers, generally under the influence of charged 
condensing agents such as multivalent ions [3, Q , short 
polyamines or charged proteins [1]. If the polyelec- 
trolytes are stiff chains, such as filamentous actin (F- 
actin) P, 0] or the fd and M13 virus 0, they aggre- 
gate into dense, hexagonally coordinated bundles. The 
counterintuitive nature of the effective attraction implied 
by bundle formation, in combination with the biologi- 
cal importance of the phenomenon, have made it the 
topic of numerous studies over the past decades (see 
Refs. d @, 0, i, i, [3 and references therein). There 
is now general agreement that strongelectrostatic cor- 
relations are a crucial condition in accordance 
with the experimental observation that most condensing 
agents carry a multivalent charge. However, other factors 
that affect the tendency of a system to exhibit condensa- 
tion are much less well established, mainly because these 
factors are difficult to disentangle in both experiments 
and theory. 

Even for a simple ionic condensing agent, its interac- 
tion strength with the polyelectrolyte is affected not only 
by its valency, but also by its size and by the sur- 
face charge density of the polyelectrolyte. In addition, 
the ionic concentration plays an important role, as it de- 
termines the stabilizing osmotic pressure exerted on the 
bundle by the surrounding solution 0, 0] and also con- 
trols the strength of entropic effects such as counterion 
release and depletion interactions Q. Thus, even a sys- 
tematic variation of the ion size can yield results that 
are difficult to interpret, as it alters both the binding of 
the ion to the polyelectrolyte and the osmotic pressure 
of the solution. Likewise, variation of the polyelectrolyte 
charge not only affects ionic binding, but also the direct 
electrostatic repulsion between polyelectrolytes. Accord- 
ingly, only an integrated approach can resolve the true 
origin of observed trends in the aggregation behavior. 



There is a remarkable dearth of such studies. Exper- 
imentally, there are technical limitations. For example, 
examination of the effect of counterion size is hampered 
by the uncertainty in measuring hydrated ion sizes [isj 
and variation of counterion valency often entails the si- 
multaneous variation of other ionic properties. Further- 
more, it is difficult to measure the ionic concentration 
within the aggregate and hence to assess the osmotic 
effects arising from a concentration imbalance with the 
bulk solution. This osmotic stress is also often ignored 
in computational and theoretical studies. Inspired by 
two recent experimental studies [13, [3]; in this Letter 
we aim to obtain a more complete understanding of how 
bundle formation is affected by (i) the size of the con- 
densing agent and (ii) the surface charge density of the 
polyelectrolyte. In Ref. the M13 virus was bundled 
using diamine molecules. In solution the diamines form 
divalent cations with a length that can be systematically 
varied. It was found that only the shortest diamines can 
induce bundle formation. This implies that increasing 
the diamine size decreases bundle stability, although a 
detailed analysis of the role of diamine size could not 
be obtained. Furthermore, it was found that increasing 
the M13 surface charge density a also destabilizes the 
bundle. In contrast, in Ref. [l3| this virus was bundled 
using alkali earth metal ions and an increased stability 
was observed upon increase of the surface charge den- 
sity It is our purpose to clarify these experimental 
findings through computer simulations. Indeed, earlier 
simulations [l3| confirmed part of the observations, but 
did not explain the surface-charge dependence observed 
in Ref. lj| and did not address effects arising from in- 
ternal degrees of freedom of the condensing agent. We 
confirm that an increase of the diamine size destabilizes 
the bundle, but also demonstrate how stability is recov- 
ered for even longer diamines through changes in their 
spatial arrangement. Furthermore, we reconcile the con- 
tradictory experimental findings for the effect of surface 
charge density. As the underlying mechanisms are highly 
generic, our findings are relevant for broad classes of sys- 
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FIG. 1: Net osmotic pressure as a function of rod-rod separa- 
tion L, for bulk diamine concentrations of 5, 30, and 60 mM, 
at a polyelectrolyte surface charge density cr — 0.633 e/nm^. 
The diamine charge separation is (a) & = 2.5 A, and (b) & = 
3.75 A. Error bars are comparable to the symbol size. 



terns that display electrostatically induced aggregation. 

We employ Monte Carlo simulations of a model based 
upon the experimental systems [l^, 14 1. M13 is mod- 
eled as an infinitely long cylindrical rod, with monova- 
lent charges placed on a rectangular grid wrapped around 
the rod at a radial distance of 28 A. The surface charge 
density is controlled via the lattice parameters of the 
grid. Soft repulsive l/r^'^ interactions between the cylin- 
der and other particles bring the effective cylinder ra- 
dius to i? = 30.5 A. The diamines are modeled as rigid 
straight molecules composed of soft repulsive beads of ef- 
fective diameter d = 5 A. The terminal beads, both car- 
rying a monovalent charge, have a center-to-center sepa- 
ration (5 and are connected by — 1 regularly spaced 
uncharged beads. The coions are monovalent beads of 
the same diameter. Water is represented as a homoge- 
neous dielectric medium (e = 80) and the temperature 
is set to T = 298 K. Electrostatic interactions are cal- 
culated via Ewald summation. The polyelectrolytes are 
placed in a periodic cell, forming an infinite hexagonal 
array with center-to-center rod-rod separation L. To re- 
duce finite-size effects a 2 x 2 array of rods is used. Rod 
degrees of freedom are ignored. As in Refs. 0, [l^], we 
employ the grand-canonical ensemble to ensure that the 
bundle is in chemical equilibrium with a bulk solution of 
diamine salt (diamines plus coions). Thus, in addition 
to a fixed number (300-500) of diamines that balance 
the polyelectrolyte charge, the bundle contains a fluc- 
tuating amount of diamine salt. We perform separate 
simulations of the bulk solution to establish its osmotic 
pressure Ilbuik as a function of chemical potential. Me- 
chanical equilibrium is then obtained if the net osmotic 
pressure 11 = (Hbundio — Hbuik) vanishes. 

Figure [11^ a) shows the net osmotic pressure as a func- 
tion of rod-rod separation, for diamines with length 
(5 = 2.5 A, at three bulk diamine concentrations. Nega- 
tive n corresponds to bundle contraction; a zero-crossing 



FIG. 2: Effect of diamine charge separation on bundle sta- 
bility for polyelectrolytes with cr = 0.633 e/nm^. (a) Stable 
rod-rod separation at fixed 60 mM diamine concentration, (b) 
Bulk diamine concentration C* at which bundling first occurs, 
(c) Separate contributions to the net osmotic pressure, for 
L = 80, 85, and 90 A, at fixed 60 mM diamine concentration. 
The dashed curves give the negative of the electrostatic con- 
tribution (i.e., electrostatic attraction), and the solid curves 
give the short-range contribution to the pressure. 



at small L indicates a free-energy minimum and yields 
the rod separation of a stable bundle. At low diamine 
concentrations 11 is always positive and no condensation 
takes place, but as the concentration exceeds a threshold 
value the pressure crosses zero and bundles form. This 
agrees with experimental results on M13 and fd [isl. IT^. 
and highlights the importance of the excess bulk diamine 
concentration. Experimentally it is found that the 
bundles become unstable when d increases. This is also 
reproduced by our simulations: for S — 3.75 A [Fig.[TI^b)] 
the osmotic pressure curves are shifted outwards and up- 
wards, reflecting a decrease in stability. To study this 
in detail, we vary the diamine length over a much wider 
range, 2.5A < S < 22. 5A. Surprisingly, two distinct 
regimes emerge. The trend observed experimentally and 
pictured in Fig. [1] continues up to (5 « 7.5 A, but for larger 
diamine lengths — not studied in Ref. |U] — this trend re- 
verses and bundle stability increases with d. 

To quantify and understand this behavior we refer to 
Fig. [21 The stable rod separation L* (at fixed, suffi- 
ciently high diamine concentration) is shown vs. diamine 
length in Fig. HJa), and the threshold diamine concen- 
tration C* required for condensation [l^ in Fig. [U^b). 
In both panels, a vertical line separates the two regimes. 
For smaller diamines, S < 7.5 A, the bundle swells lin- 
early with S. This is accompanied by an increase in C*, 
i.e., a larger bulk osmotic pressure is needed to maintain 
bundle stability. To understand this decrease in stability 
it is instructive to consider how the partitioning of di- 
amines between the bundle and the bulk changes when i5 
increases Q • On the one hand, larger diamines are more 
likely to be excluded from the bundle, thus enhancing 
the stabilizing effect of the bulk solution. On the other 
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FIG. 3: Fluid structure within stable bundles, projected onto 
a plane perpendicular to the bundle axis, for three different 
diamine lengths (top to bottom: 5 = 2.5, 7.5, and 20 A). 
Left: diamine concentration; right: diamine orientation with 
respect to polyelectrolyte axis. For discussion see the text. 



hand, swelling allows more diamine salt to enter the bun- 
dle, decreasing the concentration difference. We find that 
the latter effect dominates: as S increases from 2.5 A to 
7.5 A the diamine salt concentration of in the bundle al- 
most doubles. Accordingly, the bulk solution becomes 
less effective at holding the bundle together. 

What drives the swelling of the bundle? Wong et 



al. [1J| propose that longer diamines behave effectively 
as two monovalent ions rather than a single divalent unit, 
with the consequent loss of electrostatic attraction lead- 
ing to dissolution of the bundle. We examine this idea in 
Fig-IH^c), which shows how the separate contributions to 
the net osmotic pressure vary with diamine size. Three 
different rod separations are plotted, spanning the range 
of Fig. [SJa). The solid curves (open symbols) show the 
short-range contribution to the pressure, which includes 
the kinetic and excluded- volume terms, and the dashed 
curves (filled symbols) show the negative of the electro- 
static contribution. Thus, the intersection of two curves 
corresponds to 11 — 0, i.e., a stable bundle for the given 
6, and in the shaded regions 11 is negative. The top pair 
of curves (L = 80 A) is typical for the regime <5<7.5 A. 
Rather than the suggested loss of electrostatic attrac- 
tion we see that the electrostatic contribution is rel- 
atively constant, and it is instead the short-range repul- 
sion that changes, rising sharply with i5 and causing the 
bundle to swell and lose stability. 

We now turn to the regime of large 5. Near S = 7.5 A 
there is a weak change in slope in Fig. [D^a), indicat- 



ing a slower swelling of the bundle. More importantly, 
the threshold diamine concentration starts to decrease 
[Fig. [21b)]. This reversal of C* is particularly remarkable 
in view of the continued swelling of the bundle, as it im- 
plies that bundle stability increases despite the weaken- 
ing of the salt imbalance. Indeed, Fig. [5fc) shows that at 
large 5 the electrostatic attraction no longer remains con- 
stant, but instead rises. Important insights into this phe- 
nomenon can be obtained from the distribution and ori- 
entation of diamines in the bundle, as shown in Fig. [3] for 
fixed bulk concentration. The three rows correspond to 
S — 2.5, 7.5, and 20 A, each shown at the corresponding 
stable rod-rod separation. The left-hand column shows 
the spatial distribution of diamines in a cross-section of 
the bundle, with darker shading corresponding to higher 
concentration. The right-hand column shows how the di- 
amines are oriented with respect to the rods, with black 
indicating a tendency to align parallel to the rods (out- 
of-plane), grey indicating perpendicular alignment (in- 
plane), and white indicating isotropic orientation. For 
small diamines (S = 2.5 A) the bundle is quite compact, 
and the diamines are concentrated mostly in the two-fold 
bridging sites between pairs of rods. The orientation of 
the diamines is essentially isotropic. At S = 7.5 A — 
where the bundle is swollen and has a low stability — the 
diamines are much less concentrated in the bridging sites 
but rather form a condensed layer around each polyelec- 
trolyte. The diamines remain mostly isotropic, but have 
a slight tendency to align parallel to the rods in the con- 
densed layer and perpendicular to the rods elsewhere. As 
the diamine length increases further {S = 20 A), however, 
a more complex structure appears, where the diamines 
are once again strongly concentrated in the bridging sites 
but also exhibit a high degree of alignment. Close to the 
rods the diamines orient themselves in parallel, and in 
the bridging sites they orient themselves perpendicular to 
the rods. We have explicitly verified that in the bridging 
sites the diamines are oriented such that both monovalent 
ends lie within the condensed layer of neighboring rods. 
Thus, longer diamines act as "linkers" between the rods, 
which accounts for the rise in electrostatic attraction and 
the resurgence of bundle stability. 

Lastly, we consider the relationship between polyelec- 
trolyte surface charge density a and bundle stability. In- 
creasing cr enhances the direct rod repulsion as well as 
the coupling between the rod and the condensing agent. 
Theoretical work [l^, 111 predicts that the latter effect 
dominates: an effective attraction occurs when both the 
coupling parameter S = q^i^/n and the Manning pa- 
rameter ^ = R/fj, are sufficiently large, where q = 2 is 
the diamine valency, £b = e'^ / {AireeokBT) the Bjerrum 
length, and /i — {2Trq£-BO')^^ the Gouy-Chapman length. 
Both S and £, increase with ct, so that the net electrostatic 
attraction is enhanced at larger surface charge densities. 
There is another consideration, however: more highly 
charged rods require a higher concentration of diamines 
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FIG. 4: Stable rod separation L* as a function of polyelec- 
trolyte surface charge density a, for 5 = 10 A at 60 mM bulk 
diamine concentration. The curve demonstrates that increas- 
ing a can either increase or decrease bundle stability. 



to maintain electroneutrality of the bundle, leading to an 
increase in the repulsive short-range contribution to the 
pressure. A priori, it is not clear which of these trends — 
enhanced electrostatic attractions, or larger excluded- 
volume repulsions — dominates, and there is experimen- 
tal evidence supporting both scenarios. With diamines 
as the condensing agent, an increase in surface charge 
density from 0.303 e/nm^ to 0.343 e/nm^ dissolves the 
bundle, and a large increase in diamine concentration is 
needed to regain bundle stability [31 ■ But for certain 
divalent cations an increase in a from 0.343 e/nm^ to 
17 1 can enhance bundle stability [l3|. Thus, 



0.457 e/nm^ 

the experimentally observed trend appears to depend on 
the condensing agent. 

We shed some light on this situation by varying a in 
our simulations from 0.326 e/nm^ to 1.740 e/nm^, at con- 
stant diamine length 5 = 10 A. For this S the degree of 
alignment within the bundle is relatively low and parti- 
tioning of diamines between bundle and bulk is quite im- 
portant, so variation of L* is a strong indication of bundle 
stability. At a fixed diamine concentration of 60 mM, L* 
indeed shows a strikingly nonmonotonic variation with 
a (Fig. [4|) . For small surface charge density the bundle 
stability increases with a, as expected from the enhanced 
electrostatic coupling, yet for cr > 0.7 e/nm^ this trend 
is reversed, owing to the large number of neutralizing di- 
amines within the bundle. This reversal suggests that 
by varying a over a sufhciently wide range, it may be 
possible to observe in a single set of experiments, with a 
single condensing agent, both of the trends which have 
previously been observed only separately. 

It is worthwhile to comment on the five dimension- 
less parameters that characterize the system [l8|. The 
coupling parameter S and the Manning parameter ^ are 
convenient choices for two of these. Two further pa- 
rameters account for the role of rod-rod separation and 
bulk diamine concentration. To account for diamine size. 



Ref. [1J| proposed i(j = ^/6, with the suggestion that 
ijj > 1 is necessary for bundle formation. However, our 
demonstration that bundle stability is nonmonotonic in 
6 makes this choice of ip problematic. Moreover, the 
criterion fi/S > 1 predicts a dependence on tempera- 
ture, valency, and dielectric constant that contradicts 
well-established trends for bundle stability. Instead, we 
suggest the parameter S/ag, where is the charge sep- 
aration on the polyelectrolyte surface. This choice was 
considered in Ref. [IJ] , but disregarded because it did not 
vary significantly in the experiments. Here, we have var- 
ied 6 /as over a much larger range. Interestingly, if this 
is done through variation of a (Fig. |4]), the most stable 
condition indeed corresponds to 6/as = 0{1), i.e., the di- 
amine charge separation matches the charge separation 
on the rod surface. However, we also note that this pa- 
rameter cannot capture the nonmonotonicity observed as 
a function of S [Fig. [2Kb)], as it docs not account for the 
alignment effects occurring for long diamine molecules. 

In summary, we have clarified several mechanisms gov- 
erning the role of condensing agents in polyelectrolyte 
bundling. We have also presented a unified picture for the 
dependence on surface charge density, combining seem- 
ingly conflicting experimental observations. 

This material is based upon work supported by the 
National Science Foundation under CAREER Grant No. 
DMR-0346914 and Grant No. CTS-0120978 via the Wa- 
terCAMPWS Science and Technology Center. 
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In Ref. [l3|, variation of the surface charge was realized 
by exchanging M13 for the structurally identical fd virus. 
For each 5 the threshold concentration is obtained 



5 



through interpolation of simulation results for concen- 
trations bracketing C* , using histogram reweighting. 
[17] These surface charge densities are calculated based upon 
the value reported in Ref. [ij and the difference between 



M13 and fd. 

[18] We assume that the diameter of the diamine beads can 
be absorbed in a redefinition of the rod radius [11]. 



